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o ■ 

1 In order to study the influence of quenched disorder on second-order phase transitions, high- 

' temperature series expansions of the susceptibihty and the free energy are obtained for the 

' quenched bond-diluted Ising model in d = 3-5 dimensions. They are analysed using different 

Ch ' extrapolation methods tailored to the expected singularity behaviours. In d = 4 and 5 dimensions 

^ , we confirm that the critical behaviour is governed by the pure fixed point up to dilutions near the 

geometric bond percolation threshold. The existence and form of logarithmic corrections for the 
^•f^ ' pure Ising model in d = 4 is confirmed and our results for the critical behaviour of the diluted 

system are in agreement with the type of singularity predicted by renormalization group consider- 
ations. In three dimensions we find large crossover effects between the pure Ising, percolation and 
random fixed point. We estimate the critical exponent of the susceptibility to be 7 = 1.305(5) at 
, the random fixed point. 
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PACS numbers: 

05.50.+q Lattice theory and statistics (Ising, Potts, etc.) 
64.60.Fr Equilibrium properties near critical points, critical exponents 
75.10.Hk Classical spin models 



, 75.10.Nr Spin-glass and other random models 
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^3 '. I. INTRODUCTION 

Since many years random Ising models have served as paradigmatic systems in which the influence of quenched 
disorder may be studied through different techniques. In t he present work we rep ort results obtained by high- 
temperature series expansions. Systematic series expansions (Do mb and GreenL Il974) for statistical physics models 
' defined on a lattice provide an useful complement to field-thcorctical renormalization group studies and large-scale 
I numerical Monte Carlo simulations. This is in particular true when studying phase transitions and critical phenomena 
of quenched, disordered systems, 
^f-^ Series expansions techniques treat the quenched disorder average exactly and the infinite-volume limit is implicitly 

\^ , implied. Therefore one can obtain exact results up to a certain order in the inverse temperature for many quantities. 

' Moreover, one can keep the disorder strength p as well as the dimension d as symbolic parameters and therefore 
^^O , analyse large regions of the parameter space of disordered systems. The critical part of the series expansion approach 
' lies in the extrapolation techniques which are used in order to obtain information on the phase transition behaviour 
from the finite number of known coefficients. While for pure systems this usually works quite well, one can question 
the use of these extrapolation techniques in disordered systems, where the singularity structure of the free energy o r 
susceptibility may be very complicated, involving Griffiths-type singularities or logarithmic corrections llCardvill999l) . 
. Yet our work indicates that at least for the model we consider here, the extrapolation techniques are of comparable 
(~| ' quality as for the case of the pure (no disorder) Ising model. 

O ■ We consider in this paper the Ising model on a hypercubic d-dimensional lattice Z"* with bond dilution as a 
O ' realization of quenched and uncorrelated disorder. The pure model has a second-order phase transition for d > 2, 
and the upper critical dimension, where mean-field beh aviour sets in, is c?„ = 4. The influence of quenched disorder 
can be estimated by the Harris criterion l)Harrisl Il974j) : For d > = 4, the disorder is expected to change only 
non-universal quantities such as the transition temperature T^. At the upper critical dimension d = 4, logarithmic or 
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even more subtle correction terms should appear, and in d = 3, the phase transition of the disordered system should 
be governed by a new "random" fixed point since for the pure model the critical exponent a of the specific heat is 
positive and disorder should hence be a relevant perturbation. Our high-temperature series analyses presented in this 
paper affirm this picture. 

The rest of the paper is organized as follows: In Sect. |n] we first briefly recall the model and some of its properties. 
Section ITTll is devoted to a description of the methods used by us for generating the series expansions, and Sect. II Vl 
first starts out with some remarks on the analysis techniques used. The main results are presented in the foUowings 
subsections, where we discuss our results for the random-bond Ising model in five fSect. HV.Ap . four fSect. HV.B|) and 
three fSect. IrTU)! dimensions. Finally, Sect. IVl contains our conclusions. 



II. THE MODEL 

The ferromagnetic disordered Ising model on hypercubic lattices Z"* is defined by the partition function 

{s.} V to) / 

where f3 — l/ksT is the inverse temperature, Jy are quenched (non-negative) nearest-neighbour coupling constants, 
and the spins Si can take on the two different values ±1 . In our series expansion the combination Vij = tanh(/3Jjj) 
will be the relevant expansion parameter. Quenched disorder averages [. . ■]p{j) such as the free energy 

-/3F= [lnZ]p(j) (2) 

are taken over an uncorrelated bimodal distribution of the form 

F(J,,) = il^p)S{J,, - Jo) +pSiJ^j), (3) 

which corresponds to bond dilution: With probability p, bonds are effectively absent from the lattice, so that p — 
represents the pure system. The expansion parameter for averaged quantities is usually taken as v = tanh(/3Jo). 



A. Thermal phase transitions 

The pure {p = 0) model has for d > 2 a second-order phase transition from the disordered high-temperature phase 
to a ferromagnetic low-temperature phase. Bond dilution decreases Tc- Approaching the geometric bond percolation 
threshold p — pc the critical temperature decreases to zero (vc ^1) since above this value only finite clusters of bonds 
are present and therefore no global ferromagnetic order is possible. Figure Ogives an overview of the phase diagrams 
as calculated in this work from high-temperature series of the susceptibility. The influence of quenched disorder on 
the thermal fixed point is different above and below th e upper critical dimension du — 4. For d > 4 the pure fixed 
point is, according to the Harris criterion ijHarrij. fl974 t). stable against the influence of disorder and one expects that 
the renormalization group (RG) flow goes along the critical line from the unstable percolation fixed point straight to 
the stable pure fixed point. The critical behaviour along this line of phase transitions is governed by the pure fixed 
point. For d = 3 the pure fixed point is unstable and one expects a new random fixed point governing the critical 
behaviour in the region < p < Pc, attracting RG fiow from both sides (pure and percolation fixed point). More 
specifically, for d = 4 random disorder changes the form of the logarithmic corrections to the Gaussian fixed point 
l|AharonvLll976|) . whereas for d = b Gaussian behaviour is expected for p = and p > as well. 



B. Crossover to the percolation point 

The cross over behaviour of the diluted system near the percolation threshold of the underlying lattice is well- 
understood ijAizenman et all Il987t IConigliol Il98ll: IStauffer and Aharo nv. 1992) and independent of the dimension 
of the lat t ice. T he percolation critical point p = Pc, Tc = {vc = 1) is always unsta ble against thermal fluctu ations 
l|Coniglio[ Il98l[) and plays for our system the role of the strong disorder fixed point ijigloi and Monthu4l2nnRl) . 
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FIG. 1 Critical coupling Vc as function of the dilution p for d = 3, 4, and 5 dimensions. The arrows on the top axis indicate 
the percolation thresholds. 



The crossover exponent </> is exactly 1 and critical exponents approach their percolation values with a scaling 
behaviour of the form 

/ -2Jo/T\ 

X- (4) 
- be (^xi^), (5) 

where is a crossover-scaling function, and the critical line near v — 1, p = pc has the form 

1 — f , . , . 

—— (x{pc-p). 6) 

1 + V 



III. STAR-GRAPH EXPANSION 

The star-graph expansion method has been reviewed many times in the literature l|Singh and Chakravartvi ^^3) 
as well as in our previous papers ( Hellmund and Jankc, 2003, 2005). We will only present some algorithmic details of 
our approach and make some remarks comparing different techniques for series generation. 

"Star graph" is an old-fashioned name for 1-vertex-irreducible graphs. For observables F which allow a star-graph 
expansion (like the free energy and various inverse susceptibilities) this technique delivers a representation of F on an 
infinite (e.g. hypercubic) lattice, 

F{Z'^) ^Y,E{G;Z'^)Wf{G), (7) 

G 

as sum over star graphs (1-vertex irreducible graphs) only. Here, E{G; Z'') denotes the weak embedding number of the 
graph G in the given lattice structure (^.Martin , 1974) and Wf{G) can be recursively calculated from the knowledge 
of F{G) and F(G") for all star subgraphs G' C G. 

A property of eq. Q) important for us is that the quenched disorder average can be calculated on the level of single 
graphs as long as the disorder is uncorrelated. This is related to the fact that the method considers and counts weak 
embeddings of the graph into the lattice, i.e., embeddings where different bonds of a graph are always mapped to 
different bonds of the lattice. This is in contrast to the free embeddings used in linked-cluster expansions. Here, 
different vertices and bonds of the graph may correspond to the same vertex or bond of the lattice. Therefore in order 
to calculate the disorder average one has to classify the free embeddings of a graph according to their correspondence 
to embeddings of a collapsed graph (w here some edges of G are identified) . This is essentially the algorithm proposed 
in ijMever-Ortmanns and R,eis3 . Il999l) where the collapsed graphs are called "multiple-line graphs". This method has 
not yet been used for actual calculations of disordered systems, presumably due to its combinatorial complexity. 



TABLE I Number of star graphs with N > 8 links and non-vanishing embedding numbers on Z'*. For A*' = 1, 4, 6, and 7 only 
a single star graph exists. 
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For other series generation m ethods as, e.g., the finite lattice method which holds the world record fo r the pure 3D 
Ising model l|Arisue et gi.Ll2nO^ or methods using Schwinger-Dyson equations l|But,era and ComiLll99ij) . no practical 
generalization to disordered systems is known. 

For the generation of graphs we employed the nauty package by McKay ijMcKavL Il98l|) which makes very fast 
isomorphism tests by calculating a canonical representation of the automorphism group of the graphs. By this means, 
we classified all star graphs up to order 21 that can be embedded in hypercubic lattices, see Table |l| For each of 
these graphs we calculated their weak embedding numbers for d-dimensional hypercubic lattices (up to order 17 for 
arbitrary d, order 19 for dimensions d < 5 and order 2 1 for d < 3). For this embedding count we implemented a 
refined version of the backtracing algorithm by Martin l)MartinL Il974jl . making use of a couple of simplifications for 
bipartite hypercubic lattices Z'^. After extensive tests to find the optimal algorithm for the innermost loop, the test 
for collisions in the embedding, we ended up using optimized hash tables. 

In order to calculate the contribution W{G) of a graph G of order N to the susceptibility series, one needs the 
partition sum Z{{Jij}\G) and the matrix of all spin-spin correlation functions Mki{{Jij}\G) — Tr <5sfc,s,e~'^^^'^'^^ as 
polynomials in the N variables Jy . This is achieved by using the cluster representation 

Z 0, ^2^(C) Yl (8) 

« E2^(-) n (9) 

where the sum in eq. ((HJ goes over all clusters C C G which have only vertices with an even number of bonds and in 
eq. © the sum is restricted to all clusters Cki C G in which the vertices k and I have an odd number of bonds and all 
other have an even number of bonds. The exponent c(C) counts the number of connected components of G . A graph 
with N bonds gives rise to 2^ clusters: every bond may be present or absent. The clusters can therefore be numerated 
by A^-bit integers. These 2^ integers can be sorted in such a way that two consecutive numbers d iffer by exactly oii e 
bit (corresponding to the addition or deletion of one bond) , an algorithm known as "Gray code" ijPress ieial|,|l99l. 
This gives a simple algorithm for the calculation of the cluster sums in eqs. ((S)) and (|5J: 

- Start with the full graph C = G, store the coordination number modulo 2 ("coordination bit") for every vertex 
and count the total number r of odd vertices. Execute steps 2 and 3 from below for this cluster configuration. 

- Iterate over all other 2^ — 1 Gray codes G: 

1. Calculate the next Gray code (cluster configuration). Compared to the previous code, exactly one bond 
was added or deleted. Invert the coordination bit of the two involved vertices and calculate the change in 
r. 

2. If r = then add Jlnec ^'^ product has a factor Wn — 2vij / (1 — Vij) for every non-zero bit n in 
G. 

3. If r = 2 then add Jlnec where k, I are the two odd vertices. 

The prefactor 2'^''-^^ is taken into account by monitoring the change in the number of connected components c(G) 
in each iteration step. 

The further steps for the calculation of the susceptibility series are: 

• Inversion of the Z polynomial as a series in the {vij} up to the desired order. 

• Averaging over quenched disorder, 

resulting in a matrix of polynomials in p and v. 



• Inversion of the matrix Nki and subgraph subtraction, 
WAG)^T.kAN-^)ki-i:G'^GWAG')- 

• Collecting the results from all graphs, 

All calculations are done in arbitrary-precision integer arithmetic using the open source library GMP. For the poly- 
nomial arithmetic we dev eloped our ow n optimized C++ template library using a degree-sparse representation of 
polynomials as linked lists iKnuthLll998l) . The calculations took around two CPU years on an Opteron Linux Cluster. 



IV. SERIES ANALYSES 

The estimate of critical parameters from a high-temperature series involves extrapolation from a finite number of 
exactly known coefficients to the asymptotic form of the function. Many s uch extrapo l ation techniques have been 
developed and tested for different series and are comprehensively reviewed in l|Giitt mannl . 1 1 989() . These extrapolation 
techniques are not rigorous. They make some assumptions about the expected form of the singularity at the critical 
temperature. Field theoretic techniques like the e-expansion have similar problems. 

Usually, error estimates rely on the scatter of the results of extrapolations with different parameters (like [N/M] 
Fade approximants for different values of N and M). This may seriously underestimate systematic errors coming 
from wrong assumptions about the structure of the singularity. 

In order to get a reliable picture, we will take into account several criteria, such as 

• convergence of the analysis, 

• scatter of different approximants, 

• number of defective approximants, 

• agreement between different extrapolation methods. 

The basic meth ods we use are DLog-Pade approximation and inhomogeneous differential approximants (IDA) 
l)Guttmannl . ll989(l . In order to analyze confluent nonanalytic and logarithmic corrections, these methods are applied 
to suitably transformed forms of the series. T he parameters of th ese transformations are flne-tuned according to the 
criteria listed above, a technique pioneered in l|Adler et fl/J . ll99]]) . 

We obtain from the star-graph technique series in two variables v and p. Since the algorithm involves inversion of 
polynomials (and of matrices of polynomials) one has to find a consistent truncation criterion for the resulting series. 
The crucial observations are that a) the coefficient matrix amn of the resulting series is triangular (cf. App. i.e., 
the series is of the form 

oo n 

X{p, v)^^Y^ arnnq'^v" where g = 1 - p, (10) 

n— m— 

and b) the contribution of a graph of order N starts with terms of order q^v^ . This allows different consistent 
truncation schemes. In scheme A we calculate series which are correct up to a given order N in v, i.e., we calculate 
all non-vanishing coefficients a„i„ with n < N. In scheme B, we calculate all non-vanishing coefficients amn which 
satisfy m + n < 2N. 

Both schemes are useful in different regions of the parameter space. At small and medium disorder p we approach 
the critical point by extrapolating in for constant values of p. Therefore we use the series calculated according to 
truncation scheme A in the form 

N / n \ N 

Xiv\p = const.) oc ^ ^ a„„g" ) = XI (1^) 

n=0 \m=0 / n=0 

Near the percolation threshold p Pc this gives less and less satisfactory results. As eq. © and Fig. ^ show, 
p = const. -lines are reaching the critical line at a smaller and smaller angle. In this region, it is better to extrapolate 
along the locus w = q/v = const. This is achieved by taking the series coefiicients calculated with scheme B: 

2N 

X{v\w = const.) cx a„„u;"i;"+" ^ ^ d„w". (12) 

m+n<2N n=0 



The analysis of scheme B series brings much better resuhs than scheme A in the highly diluted region p ^ Pc- But 
the quality of the analysis is still lower (e.g., the scatter of different Pade approximants is larger) than in the regions 
with small or medium dilution. This is not surprising since we can not construct from our data a series consistently 
truncated in p. This would be necessary in order to extrapolate along the locus v = const, and study the percolation 
transition directly. Nonetheless, a reliable estimate of the transition temperature Vc is possible up to p « pc, giving 
the phase diagram in Fig. Q 



A. Five dimensions 



The analysis of the 19th order series for the susceptibility of the 5D bond-diluted Ising model is a nice warmup 
in order to get an impression of the quality of different analysis methods. Unbiased DLog-Pade and inhomogeneous 
differential approximants (IDA) along different p = const, loci give values of 7 between 1.02 and 1.04 for a wide range 
of dilutions p = . . . 0.7 with a slightly increasing tendency. 

This effect - an overestimate of 7 in a naive analysis not takin g into account confluent c orrections - is well-known 
from the 3D pure Ising model. A good illustration is Fig. 2 in ijCampostrini et all WOO'^ : whereas the results for 
many different theories in the Ising universality class (</>'' etc.) converge to 7 ~ 1.237, even 25th order series for the 
Ising model give 7 « 1.244. The remedy is the use of analysis methods tailored for confluent nonanalytic corrections. 



f{v) ^ A{v, - v)--^ [1 + A,{v, - v)^' +...], (13) 
such as the methods from l|Adler oil Il99ll) called Ml and M2. The method Ml uses DLog-Pade approximants to 

F{v) = {v,-v)f~^f{v), (14) 
av 

which has a pole at Vc with residue Ai — 7. For a given trial value of Vc the graphs of Ai versus input 7 are plotted 
for different Pade approximants and by adjusting Vc a point of optimal convergence is searched. 

The M2 method starts with a transformation of the series in v into a series in i/ = 1 — (1 — ^)^S and then Pade 
approximants to 

G(y) = Ai(l-y)^lnF(y) (15) 

are calculated which should converge to 7 as y — > 1, i.e. v Vc- These methods are especially useful when taken as 
biased approximants with a given value of 7 or Ai as input. 

Using the M2 method we find over a large range of the disorder strength p = ... 0.7 excellent convergence of the 
Pade approximants with a nearly p- independent value of Ai. As an example. Fig. |21 shows a M2 plot at dilution 
p = 0.5 from which we read off 7 = 1.001(1), Ai = 0.51(2) and determine Vc = 0.227498. Our overall estimate for Ai 
for all disorder strengths p = . . .0.7 is Ai = 0.50(5). The results of the Ml method are compatible with this but 
have poorer convergence. 

For the pure (p = 0) 5D Ising model our critical values obtained from the 19th order series with the M2 analysis 
biased with 7=1 are Vc = 0.113 425(3), Ai = 0.50(2) As usual, the errors are estimates of the scat ter of different 
Pade approximants. This should be compared with the value from Monte Carlo (MC) simulations ijLuiiten et al\ . 
11999.1 of Pc = 0.113 915 2(4), i.e. Vc = 0.113 425 0(4), which is in perfect agreement. 



B. Four dimensions 



Four is the upper critical dimension of the Ising model. Without impurities the scaling behaviour of the susceptibility 
and specific heat is t hus expected to exhibit at the cr i tical point logarithmic correc t ions. With t = |1 — T /Tc\, the 
precise form reads as (|Ballesteros et oil Il998bt iKennaL l200llKenna and Land. Il99ll Il993^ 

X ^ t-i|lnt|V3, (16) 

C - |lni|i/3, (17) 

wher e a non-singular background term in C has been omitted. For the disordered case it was shown long ago (jAharonvl 
\l<^7^ by a RG analysis that the critical behaviour is modified to take the form 



X 
C 



t exp 



D\\nt\^/^ 



|lni|i/2 



exp 



^2D|lni|i/2 



(18) 
(19) 




where D = ^/6/53 ~ 0.336. The crossover to mean-field behaviour as well as higher loop corrections 
l|Geldart and De'B elf. '1993'^ modify this unusual exponential term by factors of the form |lnt|'^i. Notice that the 
expression in eq. 1)19(1 approaches zero for t — > such that C is no longer logarithmic divergent but stays finite in the 
presence of disorder. 



1. Pure Ising model 



Due to the logarithmic corrections in eqs. and ((17|l a special treatment of the series is needed. The analysis of 
logarithmic corrections of the general form 



f{v) ^ {v,^v)-'^\\n{v,-v)\' 
is possible by l)Adler and Privmanll 19811) calculating approximants for 



F{v) = [vc — v) ln(wc — v) 



fiv) 



fiv) 



(20) 



(21) 



where one expects for singularities of the form H20|l that lim„^„^ ^(^) = ^- In what follows we call this the M3 
method. 

The estimates for the pure 4D Ising model based on our 19th order series are: a) unbiased estimate, Fig. O 
7 = 1.019(1), (5 = 0.34(1), = 0.148 6 07(10) and b) biased to 7 = 1, Fig. |2 6 = 0.42 9, 7;, = 0.148 583(3) . This 
should be compared to the analysis in fStauffer and Adlei^, '1997*) of a 17th order series l)G!aunt et all Il979^ which 
gave Vc = 0.148 588 and to MC simulations (Bittner et a/., 2002) which found f3c = 0.149 697(2), i.e. Vc = 0.148 589(2). 
While the agreement of the estimates for Vc is excellent, this analysis also shows how difficult it is to obtain reliable 
estimates for the critical exponent S of the logarithmic correction term. 



2. Bond dilution 



An extensive MC study of the si<e-diluted 4D model was done in ("Balle steros et oil ll998lJ) . The authors found 
data compatible with the theoretical prediction of a Gaussian fixed point with logarithmic corrections, but a precise 
fit of the logarithmic corrections was not possible. 

A modification of the M3 analysis method which we refer to as M3a is able to take into account the specific form 
of eq. (CHll: For 



f{v) ^ (vc — v) exp D \ ln(w, 



|l/2 



(22) 



1.01 1.015 1.02 1.025 1.03 

1 

FIG. 3 Logarithmic exponent 5 as function of the leading exponent 7 for the pure 4D Ising model via unbiased M3 approximants 
at Vc — 0.148 607. A pronounced clustering of the different [N/AI] Fade approximants around the point 7 ~ 1.019, S ~ 0.34 is 
clearly observed. 




FIG. 4 Logarithmic exponent 5 as function of Vc for the pure 4D Ising model via M3 approximants biased to 7 = 1. From the 
clustering of the [N/AI] Fade approximants one reads off Vc ~ 0.148 583, 5 « 0.429. When assessing the clustering property 
notice that the y-scale is 10 times finer than in Fig. El 



we calculate approximants to 



F{v) ^ (vc - v)i-ln{vc - v)y^^ 



f'jv) 



(23) 



which should behave as Uni„^i,^ = D/2. We made an extensive analysis of our 19th order susceptibility series 
assuming the three forms eqs. (113(1 . 1)20(1 resp. ((16(1 and 122(1 resp. ((18(1 . It turns out that, while form p3(l is not 
acceptable, both forms of logarithmic corrections allow a fit of the series data, both with a disorder dependent exponent 
S or D, respectively. The quality of the M3a fits is, however, better than for the M3 fits. This is demonstrated in 
Figs. Eland El where equal ranges of Vc {Avc = 0.0003) and of the exponent {AS = AD = 0.1) are shown. 

We interpret this as an indication for the validity of the RG prediction eq. ((18(1 . Figure [3 shows the dilution 
dependency of our estimates for the critical parameter D. It illustrates the difhculty in differentiating between a 
behaviour of the form eq. 1(20(1 and the form eq. 1(22(1 by giving a non-zero result for D a.t p = 0. It is nevertheless 




quite impressive that we can see at all such weak correction terms to the leading singularity and even estimate the 
parameters. 

For the specific heat, one expects the weak logarithmic singularity in the pure case [eq. H17(l ] to be washed out 
by the disorder. An analysis of the free energy series [F(v)]p(^j^ — T,ai(p)v^ is usually more difficult and gives less 

satisfactory results compared to the susceptibility series. One looses two orders in v by calculating C — 0^^[^ and, 
more importantly, on bipartite lattices the series for F{v) includes only even powers of v. So it has to be considered 
as a much shorter series in u^. In four dimensions we calculated F{v,p) up to order 18, which gives a series for C(u^) 
of order 8 in w^. 

Another difficulty is that the non-singular background terms are more influential in this case. Therefore DLog-Pade 
approximants show very poor convergence and one has to use inhomogeneous differential approximants (IDA) which 
are able to take background contributions into account. By this method, polynomials Pn{v),Qm{v) and Rl{v) (of 
order M and L, respectively) are determined such that 

PN{v)bf{v) + QM{v)f{v) + Rl{v) = o(i;(™+^)), (24) 

where many different triples [N, M, L) with N + M + L < (order of /) and two different variants of the differential 
operator Z?, either Di — vd/dv or D2 — d/dv, are used. The critical point Vc is then given by the smallest positive 




real root of Pn{v) and the critical exponent by (3a/(wc)/^jv(''^c)- 

We applied IDAs tailored to power-like singularities ~ t^°'^^ to the fourth derivative d'^f{v)/dv^. The results are 
shown in Fig.|Hl In the pure case we find a result consistent with d^f{v)/dv'^ ~ {vc — indicating a logarithmic 
singularity in C ^ f {v) / dv^ . At nonzero values of the dilution parameter p the singularity of the fourth derivative 
of / is clearly weaker than {vc — indicating the absence of a divergence in the specific heat in the disordered 

case. 



C. Three dimensions 



Let us finally turn to the physically most important case of three dimensions. The quest for a determination of the 
properties of the expected random fixed p oint in the SD disordered Ising i nodel is already rather long . A comprehensive 
compilation of results can be found in llFolk et 0,11 l200(i l2003(l and l)Pelisetto and Vicari 12002^ . showing a wide 
scatter in the criti cal exponent s of diffe rent groups, presumably due to la rge crossover effects. Recent MC simulations 
l)Ballesteros et a/.l . [l998a: Berche 6^"all . l2005. 20o3 llvanevko et qiJ . l2005j) provide evidence for the random fixed point 
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FIG. 9 Critical exponent 7 as function of p for the diluted 3D model obtained from about 300 DLog-Pade and inhomogeneous 
differential approximants for each dilution value. 
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FIG. 10 Critical exponent 7 as function of p for the diluted 3D model obtained from series analyses with methods Ml and M2 
(see text for details). The data of Fig. 1^1 (Pade/IDA) are also shown for comparison. 



but also show large crossover effects due to the interference with the pure fixed point for p ^ and with the percolation 
fixed point for p ^ pc (recall the discussio n in Sect. III. Bll. 

Early attempts using series expansions l)Rapar)orti Il972alb|) already indicated that the series are slower convergin g 
and more difficult to analyse than in the pure case. A review of earlier work can be found in ('Stinchcombd','1983*). 
Series analysis in crossover situations is, in fact, extremely difficult. If the parameter p interpolates between regions 
governed by different fixed points, the exponent obtained from a finite number of te rms of a seri e s exp ansion must 
cross somehow between the two universal values, and does this usually quite slowly l) Janke et all 12005)) . The mere 
existence of a plateau in 7(p), however, is an indication that here a truly different critical behaviour is seen. 

Our results are obtained from a large number of DLog-Pade and inhomogeneous differential approximants (around 
300 for every dilution value) applied to the 21th order susceptibility series x(Pi^) compiled in App. A. The resulting 
estimates for the critical exponent 7 shown in Fig. |^ however, do not exhibit any sign for a plateau. 

Since confluent corrections are essential to understand crossover situations, we also performed a careful analysis 
using the Ml and M2 methods which both take such corrections explicitly into account, see Fig. ^| This figure 
shows results using both v ^ tanh(/3Jo) as well as the coupling /3Jq as expansion variables (denoted by, e.g., Ml/v 
and Ml/ J). They give an indication of a plateau around p = 0.3 . . . 0.4, suggesting the presence of the random fixed 





point with corresponding value of the critical exponent 7 = 1.305(5). Here, as usual in extrapolation techniques, 
the error is estimated from the scatter of different approximants and we are unable to give an estimate of the 
systematic error of the ext rapolation. Our h i gh-ternperature series estimate is at least compat ible with MC results 
for site and bond dilution l|Ballesteros et ad ll998aHBerche ei"all 12004 ICalabrese oil l2003l) which cluster quite 
sharply around 7mc = 1-34(1). Field-theoretic RG estimates jFolk et al}. '2000 ^ [ Peh^e tto and Vic ariL |2000[ ) favor 
slightly smaller expo nents of 7rg — 1-32 - 1.33, while experiments iBelange r et all 11986: Birgen eau et a.l\ . fl983t 
iMitchell et all \l9H(t\ report values between 7exp = 1-31 - 1.44; for a more detailed compilation, see e.g.. Table 1 
in Ref. ijBerche et a^J . l2005(l . Since it would be extremely demanding to further extend the series expansions in the 
disordered case, better series analysis methods for the case of crossover situations would be clearly desirable. 

V. SUMMARY 

We successfully applied the method of high-temperature series expansion to the bond-diluted Ising model in several 
dimensions. The computational effort is vast and increases faster than exponentially with the order of the expansion. 
But the extensive set of combinatorical data we generated on the way, such as the list of star graphs and their 
embedding numbers into hypercubic lattices, has a large number of potential applications, for example to models with 
other kinds of uncorrelated disorder like spin glasses. On the other hand, Monte Carlo simulations of systems with 
quenched disorder require an enormous amount of computing time because many realizations have to be simulated 
for the quenched average. For this reason it is hardly possible to scan a whole parameter range. Using the method 
of star-graph expansion we obtain this average exactly. Since the relevant parameters (degree of disorder p, spatial 
dimension d, etc.) are kept as symbolic variables, we can easily study large regions of the phase diagram. 

Table m summarizes the main results of this work. Our estimates of the critical temperature for the 4- and 5- 
dimensional model are comparable in precision (as well as consistent) with the available Monte Carlo data for the 
pure case. For the bond-diluted model in 4 and 5 dimensions we demonstrated universality. The critical exponent 
7 keeps its Gaussian value 7 = 1 up to the percolation threshold of the underlying lattice, despite the fact that the 
upper critical dimension of the percolation transition is six. In 5 dimensions we also show that the exponent Ai of 
the confluent correction is universal for a large range of the dilution parameter p. In 4 dimensions we confirm the RG 
prediction that the weak logarithmic divergence of the specific heat of the pure model disappears in the disordered 
case due to logarithmic corrections of a special kind. We confirm the presence of these special corrections also in the 
case of the susceptibility. We furthermore determine the exponent of the logarithmic correction in the pure case with 
satisfying precision. 

The case of 3 dimensions is, in a sense, the most difficult to analyze. We clearly see that the pure fixed point 
is unstable and does not describe the random system. Extrapolation methods not taking into account confiuent 
corrections give a value for 7 which changes monotonously with p. By using methods tailored to the consideration 
of confluent corrections we identify a plateau in the curve 7(p) with critical exponent 7 = 1.305(5) which we see as 
evidence for the random critical point. 



model 


Vc 


7 


confluent correction 


5D 


pure 


0.113425(3) 


1.001(1) 


power-like, Ai = 0.50(2) 


disordered 




1.001(1) 


power-like, Ai = 0.50(5) 


4D 


pure 


0.148 607(3) 


1.019(20) 


logarithmic, S = 0.34(2) 


disordered 




form eq. 11811 confirmed 


3D 


disordered 




1.305(5) 





TABLE II Summary of results. In three dimensions, longer series for the pure Ising model exist, therefore we do not quote 
here our results for Vc or 7. 
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APPENDIX A: Coefficients of the 3D susceptibility 



The following table of coefHcients gives the complete information for calculating the susceptibility of the 3D bond- 
diluted Ising model up to the order v"^^ for any dilution p. 



TABLE III Table of coefficients amn of xip, v) = 1 + 6qv + 30q v + 150q v + J2 o,mnq"^v'^ up to order n = 21, where q = 1 — p. In the table, m is the column number 
starting from 4 and n is the row number starting from 4, so x = 1 + 6gv + 30q^v^ + I50q^v^ + 726g*t;* - 24q'^v^ + 3534g^u^ . . .. 
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-240 
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1344 


8400 


-1440 


-339936 


-2295744 


-12766944 


197659950 
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1608 


10296 


51480 


26544 


-1886928 


-12680496 


-70404720 


936945798 












-24 


-192 


264 


10032 


64560 


341568 


259656 


-9915696 


-69162048 


-377522064 


4429708830 










-24 


-192 


-1080 


-1704 


60024 


427920 


2062368 


2482464 


-51644200 


-367148472 


-2014331904 


20955627110 











-192 


-1080 


-12144 


-14448 


360192 


2493600 


12550416 


17926128 


-259622976 


-1931961792 


-10550435184 


98937385374 
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-1080 


-13440 


-80928 


-132024 


1840776 


14790144 


73051512 


126567264 


-1293631728 


-9980137536 


-55050628008 


467333743110 




24 
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-8544 


-95040 


-569760 


-1214880 


9797904 


82573800 


420942768 


807789264 


-6273975792 


-51221501136 


-283516855968 


2204001965006 
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1080 


5832 


-60888 


-705216 


-3910368 


-8858136 


46862760 


457439184 


2361075624 


5069434800 


-30136593768 


-259361429784 


-1455780298776 
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36720 


-437952 


-4512600 


-25012512 


-63580104 


220823568 


2449336680 


13097561328 


30177202248 


-141380350848 


-1306684851840 


-24 





1080 


24744 


129816 


283752 


-2272584 


-28419456 


-158605280 


-422656608 


936811968 


12971851368 


71258617752 


177314558064 


-655481735280 



10398318680694 
-7403140259952 
-6514909866600 



48996301350750 
-37556614417032 230940534 
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